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Abstract. We study the dynamics of a symmetric two-level system strongly 
coupled to a broadened harmonic mode. Upon mapping the problem onto a 
spin-boson model with peaked spectral density, we show how analytic solutions 
can be obtained, at arbitrary detuning and finite temperatures, in the case of 
large Q-factors of the oscillator. In general two or more oscillation frequencies 
of the two-level particle are observed as a consequence of the entanglement with 
the oscillator. Our approximated analytical solution well agrees with numerical 
predictions obtained within the non-interacting blip approximation. 
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1. Introduction 

A prominent physical model to study dissipative and decoherence effects in 
quantum mechanics is the spin-boson model [HOIS]. Currently, we witness its revival 
since it allows a quantitative description of solid-state quantum bits (qubits) [4]. A 
more realistic description requires the inclusion of external control fields as well as of 
a detector. In the spin-boson model, the environment is characterized by a spectral 
density G(uj). If the environment is formed by a quantum detector which itself is 
damped by Ohmic fluctuations, the form of the spectral density can become non- 
trivial, as it reflects also internal resonances of the detector. An example is provided 
by a flux-qubit read out by a dc-SQUID [U O [7] whose plasma resonance at ft gives 
rise to an effective spectral density G c ff(w) for the qubit with a peak at f2 [9], 
cf. Eq. ([2]) below. Recently, the coherent coupling of a single photon mode and a 
superconducting charge qubit has also been studied [10]. Until now, the effects of 
such a structured spectral density on the decoherence properties of a qubit have been 
studied in [HJ[l2] within a perturbative approach in G e s- It was shown in [13JQ3] that 
such a perturbative scheme breaks down for strong qubit-detector coupling, and when 
the qubit and detector frequencies are comparable. Hence, in [I3l[l5l[l6] the dynamics 
was investigated by mapping the spin-boson problem onto the equivalent situation of 
a two-state system (TSS) coupled to a harmonic oscillator (HO), the latter coupled to 
an Ohmic bath with spectrum Gohm- By considering the TSS and HO as the relevant 
system, analytic solutions perturbative in Gohm were obtained. 

In this work we show how to investigate the dynamics of a spin-boson system 
with a structured environment, in the case of a strong coupling between qubit and 
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detector. We evaluate the dynamics upon starting from the (nonperturbative in G e g) 
Non- Interacting Blip Approximation (NIB A) [HE]- Analytical results, valid also at 
finite detuning, are obtained by approximating the NIBA kernels up to first order 
in the detector-bath coupling strength. The paper is organized as follows: In the 
next Section we will introduce the model. Then, in Sec. [3] we discuss the well-known 
and widely used Non-Interacting Blip Approximation and its predictions. Analytical 
results for the dynamics are derived in Sec. HI 

2. The model 

In this work we consider the spin-boson Hamiltonian describing the interaction 
of a symmetric TSS with a structured environment. It reads [TJE] 

HsB(t) = ~^Y°x + ^zfiJ2~ Xk ^l + ^) + ^2hw k b~l~b k , (!) 

k k 

where er; are Pauli matrices and HA is the tunnel splitting. Moreover, b k is the 
annihilation operator of the k~ th bath mode with frequency u>k- In the spin-boson 
model the influence of the environment is fully characterized by a so-called spectral 
function, which we assume to be of the form 

G cff H = ^X k d(u 1 -Q k ) = ( ^ 2 _ cj2)2 + (r ^ )2 - (2) 

It has a Lorentzian peak of width T at the characteristic frequency fl, and 
behaves Ohmically at low frequencies with the dimensionless coupling strength a = 
lim^^o G e ff(w)/2w. As shown in [17], such spin-boson Hamiltonian can be exactly 
mapped onto that of a TSS coupled to a single harmonic oscillator mode of frequency 
O with coupling strength g. The HO iteslf interacts with a set of harmonic oscillators 
with spectral density of the continuous bath modes being Gohm(w) = kw. The 
mapping between the two models is completed with T = 2irnfl and a = 8K,g 2 /Q 2 . 
Notice that the oscillator can e.g. represent a dc-SQUID with plasma frequency fl 
which couples inductively to a superconducting flux qubit [HIE]. The damping of the 
dc-SQUID is due to its coupling to an electromagnetic environment. 

If the damping of the harmonic oscillator is small (k <C 1), as in typical 
experiments where the dc-SQUIDs are typically underdamped, then it would seem 
more convenient to use the mapping, and consider the qubit coupled to the dc-SQUID 
as unique quantum system. However, since such a system has an infinite Hilbert space, 
the inclusion of dissipation is typically done [331 [16] upon truncation of the system's 
Hilbert space to a few relevant levels (which is the case e.g. if Ml, HA ^ ksT). 
This led in [13] to find analytical results for the resonant case A = Q within a three- 
level approximation. In the present work analytical results valid at finite detuning 
A 7^ ft are obtained by focusing on the spin-boson model ^ . The advantage of this 
approach is that the reduced density matrix has rank 2. The peculiar feature of the 
peaked spectrum (EJ) is reflected in the form of the bath correlation functions, cf. (@]) 
and l|4b|l below. 

Specifically, in the spin-boson model the environmental effects are captured in 
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( 1 — cos tot) + i sin tot 



which for the effective spectral density |(2]) reads 

Q'(t) = Xt + L (e-i T costir - l) + Ze^ T sin fir + Q' Mats (r), 
Q"(r) = 7ra — e~^ T Tra (iV sin fir + cos fir) , 



being Q = y O 2 — j- and 

27T(T 



Z 
TV 



27TQ; 

n a 
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Ffi cosh (/3fl) - cos 
7ra 1 



p - fi 2 ) sinh (J3SI) + m sin 



m sinh (/3fi) + ( ^ - fi 2 ) sin 



Ffi cosh (/30) - cos 

m V 4 / 

Q'nj ats ( T ) is a function of Matsubara frequencies and has the form 



(3) 

(4a) 
(4b) 

(5) 
(6) 

(7) 

(8) 

(9) 



with the Matsubara frequencies defined as v n 
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n. For temperatures high enough 



(ksT 3> contributions coming from the Matsubara term can be neglected, as 
done in the rest of this work. 

The qubit dynamics is described by the reduced density operator p(t) obtained 
by tracing out all environmental degrees of freedom. We investigate the population 
difference P(t) := (a z )t = tr{p(t)er z }. Such a dynamical quantity P(t) obeys the exact 
generalized master equation (GME) [2] 



P(t) 



&t'K{t-t')P(t'), 



t > 0. 



(10) 



with the kernels K{t) being a series expression in the number of tunneling transitions. 

Since Eq. I|10p involves only convolutions, it can be solved by using Laplace 
Transforms. Hence, the GME transforms as 



P(A) 



1 



A + K(X) ' 



(11) 



where the same symbols P(A) and K(X) for the Laplace transform of P(t) and K(t) 
have been used, respectively. From Eq. l[TTj) . it follows that in order to obtain P(t) it 
can be enough to solve the pole equation 



A + A" (A) = 



(12) 
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and then antitrasform back to the time space. 

Due to the intricate form of the exact kernel K(t) (or K(X)), Eqs. (fit))) (or 
i[TT]0 cannot be solved neither numerically nor analytically. We must therefore invoke 
some approximations. For the symmetric spin-boson model (TJ, the so-called Non- 
Interacting Blip Approximation (NIBA) discussed in the next Section is known to 
yield reliable results over the whole regime of parameters. 

3. Non-Interacting Blip Approximation (NIBA) 

Within the NIBA [TJ [2], of the exact series expression for K(\) only the first 
term of second order in the tunneling frequency A is retained. This approximation 
has been commonly used over the whole range of temperatures and coupling strength 
to describe the dynamics of an unbiased TSS. In the undriven case, the model is 
justified for weak damping since the neglected correlations are of second order in the 
coupling a, whereas for high temperature and/or large damping extra correlations are 
exponentially suppressed. The kernel has the very simple form 

K(t) = A 2 e-Q'W cos (Q"(t)) 7 (13) 

or in the Laplace space 

/■OO 

K(\)=A 2 dre- XT e~ Q ' {T) cos (Q"(r)), (14) 
Jo 

where the bath correlation functions Q'(t) and Q"{t) have been introduced in Eq. |(4]). 
Typical results for P(t) obtained from the numerical integration of the NIBA master 
equation for the resonant case fi = A and at finite detuning ft = 1.5A are shown 
in Figs. [JJ and respectively. In the resonant case P(t) exhibits a very pronounced 
beating pattern. The analysis of the corresponding spectrum 

/>oo 

S(oj) = 2 dt cos (ujt)P(t) (15) 
Jo 

for the parameters choice of Fig.[TJ (resonant case) clearly reveals the presence of 
two frequencies, which lie around fl± « fl ± g, where g is the coupling strength in 
the TSS+HO model, as discussed above. This is in agreement with predictions of 
a three-level system Bloch-Redfield analysis (with second-order perturbation theory 
in g) as well as with exact results obtained within the numerical real-time path- 
integral approach QUAPI [13] ■ The Fourier spectrum for the detuned case in Fig. [2] 
shows a more pronounced oscillation frequency, the relative magnitude of the two 
peaks becoming larger and larger, the higher the detuning is. As one raises the 
coupling strength a between TSS and effective environment, multiple resonances 
appear, due to the fact that higher orders in Bessel functions contribute to the 
dynamics (see discussion below in the next Section). This beating pattern clearly 
originates from the peaked nature of the environmental spectrum and it is thus absent 
for the more frequently investigated cases of unstructured environments [TJ [2], i.e. 
G(uj) oc u} s e~ u '/ ulc , s > 0. Starting point is Eq. ifTTj) and its related pole equation 
{12]). The nature of the beatings as well as an analytical approximation to P(t) are 
discussed in the following Section. 
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Figure 1. Time evolution of the population difference P(t) of a symmetric TSS 
in the resonant case U = A. The parameters are: T = 0.097, a = 410 -3 , 
(g = 0.18), T = 0.1 (all quantities are expressed in units of A). In this range of 
parameters, one clearly sees that the dynamics is dominated by two frequencies. 

peaks, 
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Figure 2. Time evolution of the population difference P(t) of a symmetric TSS 
in the case of a finite detuning f! = 1.5A. The parameters are: T = 0.145, 
a = 510 -3 , (g = 0.3), T = 0.1 (all quantities are expressed in units of A). For 
a TSS being off-resonance with the HO, one notices that the relative magnitude 
of the two peaks become larger and larger, the higher the detuning is (see the 
Fourier Transform in the inset). 
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4. Weak-Damping Approximation (WDA) for a symmetric TSS 

For a symmetric TSS, NIBA is expected to be justified in the whole regime of 
parameters. In particular we have seen, cfr. Fig.[TJ that in the regime of resonance 
and strong coupling, i.e. g 3> T, it predicts the two oscillation frequencies already 
found in Ref. [T3] within a three-level approximation. However, the analysis in Ref. 
|13| was restricted to the case A = O ^> k^,T '. In the following, we shall derive an 
analytical expression for P(t) valid for arbitrary detuning |A — 0| 7^ and low-to- 
moderate temperatures k^,T < SO. The key idea is that, since we are looking to a 
sharply peaked spectral density, i.e. k = r/27rO < 1, an expansion of the NIBA 
symmetric kernel (f!4|) up to first order in k is justified. 

Since the bath-correlation functions Q' and Q" (Eq. <|4b[) ) depend in a nontrivial 
way on k, this requires some attention. In the end we obtain 

Q'oir) Q[(t) 
> * > , * , 

Q'{t) = Y (cosOr- 1) + At cosflr + Bt + C sinnr +0[k 2 } , (16a) 

Q"{t) = WsiaClT + V (l -cos Or - -TsinQrj +0[k 2 ], (16b) 

Q'o(t) 

with the zero-order terms 



4g 2 sinh/30 4o 2 , s 

Y = — W = — (17) 

Q 2 cosh/30 - 1' o 2 ' K ' 



and first-order terms 



A = - T ^> ( 18 ) 

2 5 2 ^0 + 2sinh^O 
C " L fi3- cosh/30-1 ' (2Uj 

v^r^. (2i) 

Notice that the contribution coming from the Matsubara frequencies iJH) has been 
neglected, affecting only the short-time dynamics. We will here discuss first the simpler 
undamped case (k = 0) and later perform the weak-damping approximation on the 
NIBA kernels. 



4-1- Undamped case (k = 0) 

In this Subsection we discuss the case of a TSS coupled with an undamped HO 
initially prepared in a thermal equilibrium state. The pole equation reads now 

X P + K (X P ) =0, (22) 



with 
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Ko(X) = A 2 / dre- Ar e~ Q '^ cos (Q' '(t)), 
Jo 



(23) 
(24) 



where we denoted with X p the solution of the undamped pole equation. Notice that 
-Ko(A) has the same expression as in Eq. (fl4|) . if one replaces Q' and Q" with Q' and 
Qq, respectively. In order to evaluate Eq. (|23| analytically, we replace cos (Qq (t)) 
with 5He{exp (iQ' ' (t))} and we perform the Jacobi- Anger expansion [18] 



J Q (z) + 2 ^2 [n Jn(z) cos (ny) , 



(25) 



where J n (z) are Bessel functions of a complex argument. These expansions are valid 
only in the case in which z is independent of y, which is the case for a TSS coupled 
to an undamped HO. We also make use of Graf's Addition Theorem 



Jn+k(u)J k (v) c s fl(ka) = J n (w)lfl(nx), 



k= — oo 



where 
and 



•\Ai 2 + v 2 — 2 u v cos a 



u — v cos a = w cos x, 
wsina = uisinx. 



We finally obtain 



K (X) = A 2 e Y I dre" AT ^e 



Jo(wo) + 2 Y i n Jn(u ) cos [n(Qt - x)\ 



where (see | Appendix A [ 



iY 



u 



cosx 



Wy 2 - W 2 



Ag 2 



X = 7T + 1- 



tanx 



n 2 



,w 



sinh 



Y 



(26) 

(27) 

(28a) 
(28b) 

(29) 

(30a) 
(30b) 



After expanding the cosine which appears in the Eq. ([29| and after noticing that 
Jo(uo) and i"J„(uo) are always real, the espression for the symmetric kernel in the 
undamped case finally reads 



POO 

K (X) = A 2 e Y / dre 
Jo 



Xt 



Jo{ u o) + 2 ^(— i) n J„(iio) cos (nflr) cosh ( n 



(31) 
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It is useful to introduce here some amplitudes, in order to enhance the readability 
of the kernel. We therefore define 



A„ (c) = Ae 172 ^ (2 - 6 n ,o) (-i)« J„(u ) cosh JrJ^j , 



(32a) 



such that we can rewrite Eq. f31j) in the very compact form 



+°o poo 

K (X) = A 2 l(c) / dre- XT cos (nfir). 



The population difference Pq(X) in the undamped case becomes 
1 



PoW = 



A + K Q (X) 



1 + A 2 - 

A 2 nr=i(A 2 +« 2 o 2 ) 



n~ (* 2 + « 2 « 2 ) + e£=o ^ (c) n:r (A 2 + 



(33) 

(34) 
(35) 

(36) 



and it is clear that the pole in A = is not a physical one, since Pq(X = 0) vanishes. 
This means that the dissipation-free (n = 0) pole equation reads 



X p + K a (X p ) = 



J (A 2 + n 2tf) + ^ Al l(c) n (A 2 + nW) = 0. (37) 



n=0 



m=0 



n=0 



4-2. The weak-damping population difference P{t) 

The weak-damping kernel Kwda(X) is obtained from Eq. lfl4|) by retaining only 
terms up to first order in the linearized in k bath correlation functions Q[ and Q'{ . It 
reads 

POO 

K WDA (X) = A 2 / dre- Ar e-«iW {cos(Q '(r)) [1 - Qi(r)] - sin (Q^tW/MI- 
./o 

(38) 



The WDA kernel will be used in Eq. (fl2|) to solve the pole equation and finally obtain 
Pv/DA{t). Consistent with the previous prescription k -C 1, we can expand the 
solutions A* of the pole equation around the solutions X p of the non-interacting pole 
equation up to first order in k In other terms 

A* = X p — Kjp + iKip, (39) 

where A p satisfies the undamped pole equation <f37|l . By inserting Eq. (fl6|) and (j39"l) 
in Eq. (fSHJ, one finds the following expressions for the kernels (to first order in k) at 
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K WDA {X*) = A 2 dTe- A " T e- Q ° (T) 
Jo 

x {cos (Q{,'(r)) [1 + k 1p t - lKcpr - Q[(t)} - sin {Q' \t))Q'{{t)} + 0[k 2 }. 



According to Eq. I|40p , the pole equation lfT2|) now reads 



-«7 P + iiitp + A / dre 



-A„T . 



(40) 



(41) 



x {cos (Q '(r)) [k 1p t - \K(pT - Qi(r)] - sin (Q' \t))Q'{(t)} = , 



where we used the pole equation for the undamped case (|22|) . After isolating the real 
and the imaginary terms from the above equation [J we find 



'Kip 



IKLp 



1 - A 2 / dre- A » T e-QoM cos (Q' Q '( T )) r 
Jo 

/>oo 

= A 2 / dre- A ^ e -^M[ cos( Q^ (T)) Q; (T)+sin( g^ T)) Q' i ' (r) ] ; 
1-A 2 / dre- A " r e- Q ° (r) cos (Q^(r))r =0. 



(42a) 



(42b) 



If the term between brackets is different from zero, one easily gets <p = and, after 
some rearrangements, 

A^dre"^ e^oM [ CO s (Q '(r)) Q[(r) + sin (Qg(r))Qi'(r)] 



7 P = 



1 + ^AVA) 



(43) 



Once we have obtained the expression for the decay rates 7 P corresponding to 
each pole A p , we have all ingredients to get the population difference P(t) with the 
help of the Residue Theorem. In fact, it holds 



P(i)^£e At P(A) 



— K7p(A p )t 



lim [A — (X p — K^p)] 

A=A — K7^Ap — K7 P 



(44) 



X + K\YDA (A) ( 45 ) 



as follows from Eq. (j40|) . Notice that here we split the damping-dependent and the 
damping- indipendent contributions as A = Ao — ivy. The limit lim can 

\=\o—Ky—>\p—K'Yp 

also be rewritten as lim lim . Hence, performing first the limit over the decay 



Ao->Aj, 7^7p(^o) 



rate 7, we find 
P{t) = e^e-^ lim (A - X p ) 



A + A 2 / dre- AoT e^'^ cos (Q '(r)) 



-K 1- A" 



dre -Ao. g-Qi(r) CQS (Q»( r )) T ) ( 7p(Ao ) _ >(Ap)) 



(46) 



I Note that the Laplace Transform of an odd function of r is even in A and vice versa. In this case, 
the integrand is odd in t, thus the corresponding Laplace Transform is even in A: For pure-imaginary 
values of A, the result of the integral is real. 
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P(t) = V e V e -"-M lim (A - A p )P (Ao) , (47) 

z — - A — *Ap 

as it follows from Eq. l(34| . We stress that the dynamics in the (weakly) damped 
case is essentially determined by the correponding undamped dynamics, the damping 
being only responsible for pole-dependent exponentially decaying factors. 



4-3. Series expression for the weakly-damped symmetric kernel and decay rate 

Now we want to find a compact analytical form for the kernel Kwda(X) and the 
decay rate j p (|38| and (j43j) . respectively. To this extent, let us start from the kernel 
Kwda(X), the generalization to the decay rate being straightforward. 

As in the undamped case, in Eq. f38|) we replace cos(Qq(t)) with 
9te{exp (iQo ( T ))} an d s i n (Qo( T )) with Jm{exp (iQg (r))}. Analogously to the 
procedure followed for the undamped kernel |33|) , by using the Jacobi- Anger expansion 
l|25"Tl we obtain 

+00 „oo 

K WDA {\) = J2 d ™- AT {A 2 n{c) co S (nfir) [1 - Q[(r)] + sin (nJ2r)Qi'(r)| , 

(48) 

where uq and x are given by Eq. (|30a[) and Eq. <|30b[l . The dressed tunneling elements 
A„ (c) have been already defined in Eq. (|32a[) and 



A„ (a) = Ae Y/2 ] j(2 - 5„, ) (-i)» J n {u ) sinh (n^j . (49) 

The expression for P(X) follows from l(48| (see the discussion in the Section l4~5"l below). 
Along similar lines, the decay rate 7 P , cf. Eq. (|43)l . may also be written as 



7 P 



1 EnZ /o^dre-V Al M cos (nOr) Q' x (r) + A» sin (nfir) Q?(r) 



v + oo A2 ZA P 
Z.„=0 »<c)( A 2 +rl 2f}2)2 



(50) 



T/ie case n = 0, n = 1 

In order to obtain a more useful analytical expression for P(t), we notice that 
the lowest orders in n give the largest contribution to the sum in Eq. 1|48|1 . because 
the amplitudes A^ depend on Bessel functions J n {x), which roughly behave as x n as 
soon as the argument becomes small. Since we investigate a regime of temperatures 
generally smaller than fl, i.e. (3MI/2 > 1, and of coupling such that in general g < Q, 
then the quantity uq is smaller than one. Hence, we just restrict our analysis to the 
orders n = 0, n = 1 in Eqs. I|48p and (|50| . The undamped pole equation Eq. (J37j 
therefore becomes (we identify here A with A 0(c) for the sake of clarity) 



(xl + n 2 )(\l + Al) + x 2 p A 2 1(c) =0, 



(51) 
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n 2 



± 



\ 




(52) 



We notice that only terms quadratic in X p appear in the formal expression of the decay 
rate (cf. Eq. l|B,l|0 . Hence, it is enough to express the poles as in Eq. (|52|) , 

Given the poles in the undamped case, we can substitute each of them in the Eq. 
|50| for 7 P with sum restricted to n = 0,n= 1. We will refer to them as 7± = 7(A±), 
the explicit form of the decay rate being given in |Appendix B| 



4-5. Analytical expression for P(t) 

In order to obtain the analytical expression for P(t) in the symmetric case, let us 
start again from Eq. lj47|) . By summing up all residues contributions, we end up with 



P(t) = 



— p-K-1-t 



A 2 _ + 
A 2 _ + A 2 



cos 51_t 



\ 2 

,-K-t+t 



Q 2 



A 2 + A 2 _ 



a 2 _ + n 2 



A 2 + 



sin Q-t 



cosfi+t 



n + x 2 + -x 2 _ +t > 



(53) 



where Q± = — i\±, as follows from Eq. (]52[) . 

Notice that the expression for P(t) is invariant upon exchanging the frequencies 
Q- -► Q+. 

The analytical formula Eq. (f53|l for P(t) is compared in Fig. [3] with the exact 
numerical NIB A and the conventional weak-coupling approximation (WCA), obtained 
by a linear expansion of the bath correlation function for the coupling strength a [2]. 
There, the analytical form for the probability difference reads in the symmetric case 



P{t) = jcos At + sin At | , 



-7*>* 



(54) 



2k B T 



where 7^ = ^5(A) is the dephasing rate. Moreover, S(lo) = G e &(cj) coth ^ 
a spectral contribution which represents emission and absorption of a single photon. 
The choice of parameters in Figs. [3] and [4] is the same as in Fig. Q] and as Ref. [13], 
under the resonance condition O = A. One can notice a very good agreement between 
NIBA and the analytical WDA, whereas Eq. (|54|) completely fails in describing the 
oscillatory behaviour of P(t). In Fig. |4] the corresponding Fourier Transform of the 
probability difference is showed. There, one can see the missing oscillation frequency 
of the conventional WCA given by Eq. (|54")l and the excellent agreement between the 
numerical NIBA and our analytical solution WDA. 

Finally, in Fig. [5] we show a comparison among the WDA and the NIBA in 
presence of finite detuning |A — fl| = 0.5 for a higher coupling strength between qubit 
and HO (g = 0.3), keeping the coupling between detector and environment constant. 
Also in this case, the weak-damping approximation fully agrees with the numerical 
solution of the NIBA. In the inset one can also notice the disagreement of both models 
with the conventional WCA, characterized by a single oscillation frequency at u = A, 
see in particular Fig. [6] which shows the Fourier Transform of P(t). 
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— WDA 

NIBA 
- conventional WCA 



Figure 3. Time evolution of P(t) within the NIBA as well as the analytical 
WDA. The parameters are as in Fig.Q] namely = A, V = 0.097, a = 410~ 3 , 
(g = 0.18), T = 0.1 (in units of A). Notice the perfect agreement between the 
numerical NIBA and the analytical WDA. A perturbative approach in G e s(w), 
denoted here as "conventional WCA" (see J54J), completely fails because it does 
not even account for the two main oscillation frequencies. In the inset the short- 
time dynamics is magnified. 
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— WDA 




i NIBA 






- conventional WCA 









1.5 



0) [A] 



Figure 4. Spectral function of P(t), corresponding to the same regime as in 
the previuos case. One clearly sees that the Fourier Transforms of NIBA and 
WDA exhibit a double peak structure. In contrast, the WCA predicts a single 
broadened oscillation peak. 



5. Conclusions 



In conclusion, we discussed the dynamics of a symmetric TSS interacting with 
an effective structured environment, modelling a qubit interacting with a readout 
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t[A] 



Figure 5. Time evolution of P(t) at finite detuning within the NIBA and the 
analytical WDA. The parameters are as in Fig. [2] namely Q = 1.5A, T = 0.145, 
a = 5 10~ 3 , (g = 0.3), T = 0.1 (in units of A). Again, the agreement between the 
numerical NIBA and the analytical WDA is striking. The conventional WCA l|54ll 
(see inset) also in this expected, fails in describing the correct dynamics. 
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Figure 6. Spectral function of P{t) for the NIBA and the analytical WDA (same 
parameters as in the previous Fig.[5j. The oscillation frequencies of WDA and 
NIBA coincide, whereas in the inset one can at first glance see the appearance of 
a single peak concerning the Fourier Transform of the conventional WCA. 



dissipative detector. This case has not been so far deeply investigated within an 
analytical approach out of resonance, i.e. if the tunneling frequency A differs from 
the detection frequency f2. We approximated an exact generalized master equation 
(GME) within a novel weak-damping approximation (WDA) which, in contrast 
to "conventional" weak-coupling approaches [2], is able to correctly reproduce the 
dynamics, characterized by multiple oscillation frequencies. The WDA approach is 
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based on a first approximation of the kernel of the GME up to second order in the 
tunneling frequency, i.e. the Non-Interacting Blip Approximation (NIB A), which for 
a symmetric spin-boson model is valid over the whole range of parameters. Then, 
in order to obtain an analytical form for the dynamical population difference Pit), 
for small enough temperatures (i.e. k^,T < Ml) and small coupling strength between 
detector and environment, we approximated the NIBA kernels within a weak-damping 
approach, whose details are explained in Sec. [H The agreement of our analytical 
solution for P(t) valid at arbitrary detuning |A — fi| ^ with the numerical NIBA 
is striking. The former one is able to reproduce the two oscillation frequencies which 
are related to the tunneling and the detection frequency, respectively, as predicted 
by ab-initio numerical schemes like QUAPI [13]. Our results are of interest for the 
understanding of dephasing in qubits strongly coupled to a broadened harmonic mode 
as e.g. flux qubits [HE] or cavity QED qubits [10] . 
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Appendix A. Calculation of the parameter x 



Here we are interested in finding out the correct value for x which comes in the 
formulas after we make the expansion in Bessel functions. We refer here to the Eq. 
(|38l) . which we rewrite for clarity. 



dre -Ar e -Oj(r) {CQS (Q « (r)) {1 _ Q / (r)] _ sin (Q»( t ))Q'{(t)} . (A.l) 



K(X) 



Let us examine, to fix the ideas, the term exp{— Q' {t)} cos (Qo(r)): 
r Q '° {T) cos (Qo(r)) = e Y me{e- YcosnT e +iWsinnT } (A.2) 

^ cos Qt-\-W sin Qr] ^ ^ 

iY W 



WW 2 -Y 2 



e Y V\t{e 



It could be now convenient to interprete 

iY 



y/W 2 - Y 2 



cos SJr+ 



y/w 2 - y 2 



: sin Or 



Y 



VW 2 - Y 2 VY 2 - W 2 
-W iW 



smx 



y/W 2 - Y 2 \JY 2 - W 2 
so that the exponent can be rewritten as 

e- Q oM cos (Q«(r)) - e Y 9\c {jVW^Y* cos (Ot +x )^ 



(A.4) 

(A.5a) 
(A.5b) 

(A.6) 
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At this point we can use the Jacobi- Anger expansion |25|) to expand the exponent 
in series of Bessel functions. We finally obtain 



K(X) = A V jf "dre - ^ j 



+°° / 
Jo{uq) + 2 i) n J n (uo) co s (nOr) cosh I n—— 



[1 - Q1(t)] + 2 J2(-i) n J n (u ) sin (nfir) sinh (»^) Q'/(r)} , 



n=l 



(A.7) 

which coincides with Eq. l(48|) . once we introduce the amplitudes. Here, uq is given by 



VV 2 - r 2 = ivV 2 - vk 2 = i% — L— , (A., 

K (2 ) 



u 

sinh 



since Y" = -VKcoth(^), hence \Y\ > \W\ (note that W > and hence F < 0). One 
notices that the argument of the Bessel functions is small, whenever /3hQ/2 > 1 and 
g < f2, i.e. in the regime we are interested in. 

We would like now to obtain the exact value of x, which is easily performed. Let 
us start from Eq. (| A.5[) and let us rewrite the tangent as 

tanx = y = — itanh t, - ^ - ) = tan (— • (A-9) 

We assume x to be complex, therefore we write it as x = a + ib. In general, it 
holds 

cos (a + ib) = cos a cosh b — isinasinhfe (A. 10) 

sin(a + i6) = sin a cosh& + icos a sinhfe . (A. 11) 

Prom (|A.10[) . in order to have cos (a + ib) = -\-Yj\lY 2 — W 2 , namely a real number, 
it must be 

a = nir. (A. 12) 

From Eqs. ([AlO)) and ([All]) , we can write the tangent as 

. .,. tana + itanhfe , , ,W , . _ . 

tan (a + ib) = - : — > +itanhfe = +i— , (A. 13) 

1 — l tan a tanh b a=mr Y 

as we get by calculating the tangent from Eq. (|A,5jl . Hence, 

W 

tanh6=y. (A. 14) 

We can eventually write x as 

x = a + ib = mi + i arctanh— — nir — i—^- . (A. 15) 

In order to decide, whether to assume n = or n = 1, one must look at the cosine 
or sine: 

cos, (-D" , 1 2 = (-ir^=f = (-irA (A - 16) 

VI - tanh 2 6 / VK 2 % /y2 _ ^2 

y 2 
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+Y 



VY 2 - W 2 
or, equivalently, 



Eq.JXTTJ 

sin a; > i(— 1) 



< 



n = 1 



tanh 6 



yl - tanh 2 b 



= i(-l)' 



W 

Y 



Y 2 



= i(-l) T 



VY 2 - w 2 



-i(-l) r 



W 



AW 



VY 2 - W 2 VY 2 - W 2 



(A.17) 



(A.18) 



n = 1 . 
(A.19) 



Appendix B. Explicit form for the decay rate 7 P 

In this Appendix we wish to give the analytical result for the decay rate 7 P (A P ) 
as function of the solution X p of the undamped pole equation Eq. lj5lj) : 



7p(A p ) 



with 



and 



K 2A 2 



(A 2 + n 2 ) 2 + A 2 ft 2 



+ Q 2 Is + wAl M + tAi 1 + X p ( Af (c) 5 (A P ) + h(X p ) 



Al4 + 2B-D) + tf(2B-^), 



p = (2A- B)n 2 A 2 + (B - D)Al 

q 

t ee --^(^ 2 + 3A 2 ), 
V 

u = — 3 



2 

IT 

l2 a 2 



W ee A 2 (|+B)-tt 2 | 



2\2 



/i(A p ) 



(A 2 + » 2 ) 
A 2 + 4ft 2 



A X 2 - 4» 2 
2" A 2 + 4ft 2 



(A 2 +ft 2 ) 2 vn 3A 2 + 20ft 2 
A 2 +4ft 2 1 A 2 +4ft 2 ' 



(B.l) 

(B.2) 
(B.3) 

(B.4) 

(B.5) 

(B.6) 
(B.7) 
(B.8) 

(B.9) 
(B.10) 



As already seen, the physical poles are A 2 = — X\ 2 = A^. Correspondingly, the 
decay rates j± = 7(A±) follow according to Eq. IjB.ip . 
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